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Abstract 

We study the entanglement and Renyi entropies of two disjoint inter- 
vals in minimal models of conformal field theory. We use the conformal 
block expansion and fusion rules of twist fields to define a systematic ex- 
pansion in the elliptic parameter of the trace of the n— th power of the 
reduced density matrix. Keeping only the first few terms we obtain an 
approximate expression that is easily analytically continued to n — > 1, 
leading to an approximate formula for the entanglement entropy. These 
predictions are checked against some known exact results as well as against 
existing numerical data. 

1 Introduction 

A useful quantity to study the phenomenon of quantum entanglement in ex- 
tended systems with many degrees of freedom is the von Neumann or entangle- 
ment entropy. It is defined as foUows. If a pure quantum state j^P) (typically 
the ground state) can be subdivided into two complementary subsystems A and 
B one can construct the reduced density matrix 

PA-trs|^')(*| (1) 

by tracing over the degrees of freedom of B. The entanglement entropy Sa is 
simply the von Neumann entropy associated to pA 

Sa = -tr PA In PA ■ (2) 

The entanglement entropy has been extensively studied in low dimensional 
quantum systems as a new way to investigate the nature of quantum critical- 
ity [IHSIIH]- Several different calculations based on the conformal field theory 
(CFT) describing the universal properties of the quantum phase transitions in 
1-1-1 dimensional systems, like spin chains, have shown that the entropy grows 
logarithmically with the size £ of the subsystem A as [2HZ] 

^A-^log^ + fc, (3) 
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where c is the central charge of the CFT and fc is a non-universal constant 
related to the ultraviolet cutoff. 

As an extension of the von Neumann entropy one considers the Renyi en- 
tropy, defined as 

4"^ = ;^logtrpl . (4) 

The von Neumann entropy can be reached in the limit n — > 1. 

The quantity tr p\ plays a major role in the replica approach to entanglement 
entropy [2j[6]. This method is based on the fact that, for integral n, tr is the 
ratio Z„(A)/Z", where Zn{A) is the partition function on a n— sheeted Riemann 
surface, obtained by joining cyclically the n sheets along region A, and Z is the 
partition function of a single sheet. In one-dimensional quantum systems the 
subsystem A consists of one or more disjoint intervals A ~ AiU A2U . . . and the 
trace tr is proportional, at criticality, to the n—th power of the correlation 
function of local primary fields sitting on the end points Ui, Vi of Ai 

trp^^cx (rK,wi)r )...)" . (5) 

These primary fields have conformal weight [5] 

A„^A„^^(l-\) (6) 
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hence do not belong in general to the Kac table. They can be considered as a 
special kind of twist fields, called branch-point twist fields [THI, because they are 
naturally related to the branch points in the n-shceted Riemann surface where 
the system is defined. 

When A consists of a single interval of length £ the equation ([S]) yields 

trp:^ = c„r (7) 

where c„ is a non-universal constant. This expression can be easily analytically 
continued to any real or complex n and the limit 

- lim dn trp'X =■ Sa (8) 

gives at once Eq. ([3]). 

When the subsystem A consists of more than one interval, the analysis be- 
comes more complicated and a complete description is still lacking. In the case 
of two disjoint intervals A ^ AiU A2 = [ui, ui] U [1*2, W2] global conformal in- 
variance gives 

/ I II I \ T^in—l/n) 

trpl^clf- "^-"f^--- X M.), (9) 

V|U1 - Vi\\U2 - V2\\Ui - V2\\U2 - Vi\ J 

where x is the cross-ratio 

(Ul - Vl)(u2 - V2) 



{Ui - U2){vi - V2) 



(10) 



It was formerly assumed [B] J-n{x) = 1 identically, but it was subsequently 
realized [T^IT^ through analytical and numerical means that this choice was 
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incorrect. This observation generated in last years an intense research work 
aimed at determining analytically and/or numerically the function !Fn{x) 
[23]. In particular in [231 been shown that in any CFT Tn{x) admits a 

small X expansion and the first few terms have been evaluated. 

Our goal in this paper is to describe a general method to calculate the 
function T„{x) in minimal models using conformal blocks and fusion rules of 
twist fields. In section [3| we apply an idea of Zamolodchikov of expanding the 
conformal blocks with respect to the elliptic variable q (see next section for 
details). Since even very small g's can be related to large cc's, this expansion, 
even by just considering only the first few terms, gives a good approximation 
for Tn{x) (see Eq. ([3T|) ). that is easily analytically continued to any n. This 
leads to write an approximate expression for the contribution of J-"„ to the von 
Neumann entropy TvNix) = lim„^i 9„J>j(a;), which fits nicely to the numerical 
data of the critical Ising model. Our formula for Tn{x) also compares favorably 
with exact results at n = 2 which are described in the last section. An appendix 
describes some properties of the two functions S2(n,a) and 54(71, a) which give 
the contribution of the two spin and four spin operators to J-"„(x). 

2 Four point function of twist fields and confor- 
mal block technique 

In this section we show how one can calculate the four point function of twist 
operators on a sphere by using the structure constants and conformal blocks. 
In the first subsection we introduce the conformal block expansion of four point 
function and the Zamolodchikov's elliptic recursion relation to calculate the 
conformal blocks. In the second subsection we summarize the fusion structure 
of the twist fields for different copies of conformal field theories. In the last 
subsection, combining the results of the two subsections, we propose some good 
approximate results for the Renyi and Von Neumann entanglement entropy of 
minimal models. 

2.1 Zamolodchikov's recursion relation 

Using global conformal symmetry one can always write down any four point 
function as 

< 01(21,51)02(22,^2)03(23,^3)04(^4,^4) > = 

( ) . ( - ) ">3^l(x, ^), (11) 

VZ14Z23Z12234/ VZ14Z23212234/ 

where 

212^34 _ Z12234 



Z13224 Z13Z24 



(12) 



is the cross ratio and is the conformal weight of the operator Oi(zi, zi). 

It is well-known that iFl^{x, x) in any conformal field theory can be expanded 
with respect to the conformal blocks as [TT] 

Fll{x,x) = Y,C[2CUF[i,Ai,^,,x)F[i,Ku~^r.3:); (13) 
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where c is the central charge of the conformal field theorj0, F{c, Ai, Ai,x) is 
the conformal block, C'2 is the structure constant and the indices I indicate 
the fusion channels. For those cases that at least one of the Oj's has a null 
vector^ one can usually write a differential equation that F{c, A; , , x) satisfies. 
Although in limited cases one can find the conformal blocks explicitly [11] it is 
usually very difficult to find a compact formula for the conformal blocks. When 
none of the Oj's has a null vector one can simply write an expansion of the 
conformal block with respect to the cross ratio. However, as was already noticed 
in |llj . this expansion converges very slowly and so one needs to do numerical 
calculations to get reasonable results. A recursion relation formula was found by 
Al Zamolodchikov [25j which is more suited for numerical calculation, however 
still the convergence is slow and one can not get interesting results by just 
taking the first few terms. Al Zamolodchikov was able to find another recursive 
formula by expanding the conformal blocks with respect to the elliptic variable 
q; P51[?r| . Since even very small g's can be related to large x's, this expansion, 
even by just taking few terms, gives very good approximation of the conformal 
blocks for large a;'s. The Zamolodchikov's formula for the conformal block has 
the following form 



F(5,Az,A„x) = {16q)'^'-'^x'^{l-x)'i^{93{T))'^-*'H{c,Ai,X,q); (14) 
i?(c,A,A.,,) = 1 + ^(16,)- ^'^f " + '^^ ; (15) 

■m,n ^mn[C) 

where <5 = X]i=i ^'^'^ 

c - 1 (Bm- p-^n)^ , , 

A™„(5) - f -■ (16) 

The q is the elliptic variable and has the following relation with the cross ratio 



' K{x) " K{x) ' "~e|(r)' 

where Kix) is the complete elliptic integral of the first kind and 0i(r)'s arc the 
Jacobi elliptic functions defined as 

1 

(t(l-t)(l-xt)) 
e2(r) = ^g("+l/2)^ 03(r)-^<^"^ (19) 

i?,„„(c, Ai) has the following complicated form 



1 ' 1 

i?,„„(c. A,;) "^IIt" 



2 Ah] 

k.l 



11(^1 + ^2- ^)(-Ai + X,- ^)(A3 + A4 - ^)(A3 - A4 - ^); (20) 



^We reserve c as the central charge of one copy of minimal CFT in the next pages. 
^In other words when at least one of them is part of the Kac table. 
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where 



P = 1 ((1-5)1/2 + (25 -5)1/2), (21) 



/24 

A C — 1 , o 

A. = ^ + A?, 
and the products are taken over the following sets: 

p = — TO + 1, — 7Tl + 3, TO — 3, m — 1, 

q = — n + 1, — n + 3, 71 — 3, n — 1, 

k = — TO + 1, — 7Tl + 2, TO — 1, 771, 
I = — 77. + 1, — 77 + 2, 7T, — 1, 71. 

The prime on the symbol of the first product in (|20p means that the factors 
with (fc, I) = (0, 0) and (777, 77) must be omitted. 

To have an idea about the expansion we write the first few terms oi H{c, A; , , q) 



as 



H{£,Ai,A,,q) = l + J2hk{c,Ai,Ai){16q)^ (22) 

k=l 

ft,(e,A,.A.) = (23) 

There are some comments in order: 

1. Looking to the equation (jf 7p one can easily see that the expansion with 
respect to q is converging quite faster than the expansion with respect to 

X. 

2. The expansion ([T5|) has a singularity whenever A; is part of the Kac table. 
This means that the expansion is not useful in calculating the four point 
function of the operators in the Kac table. 



The expansion ((T4| was derived in [27] for c < 1 and c > 25. The technique 
is based on using the expansion of classical block for large central charges 
and large conformal weights. In calculating the classical block one needs 
to study the behavior of the five point function ( the original four primary 
operators and an operator with the second level null vector). The second 
level null vector is just present in conformal field theories with c < 1 and 
c > 25 . Although in deriving the classical behavior of conformal block 
one needs to use second level null vector, it is quite likely that the final 
result be correct for generic central charges by analytical continuation. 
Another important thing to notice in this direction is related to possible 
complex numbers ( with non-zero imaginary part) that can appear in (jl4p 
for 1 < c < 25. It is easy to see that Am„(c) is not real in the region 
1 < c < 25. Using (PH)) it is not difficult to show that in this region 
Rmn{c, Ai) has also imaginary part. However, one can show that in the 
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first few terms of the expansion the imaginary parts all disappear and we 
have a real number. We assume this is true for all the terms. Since in this 
work we deal with n copies of CFT's we will usually have 1 < c < 25, so 
we will assume that Zamolodchikov's formula is correct also in this region. 

4. Based on the associativity of the conformal algebra one expects crossing 
symmctr}lf| for the four point functions of primary operators in any CFT. 
In other words 

Tliil -x,l-x)= Tliix, x). (25) 

It seems that it is a very difficult task to show analytically that the 
Zamolodchikov's formula satisfies this relation. However, it is not difficult 
to show numerically that by keeping the first few terms of the formula 
the equality is approximately valid. 



2.2 Fusion rules and structure constants of twist operators 

To calculate the four point function of the twist fields by using conformal block 
technique one needs to know the fusion structure of two twist fields. The OPE 
of two twist fields can produce many different operators, in case of minimal 
models for generic n we do not know how to write all the primary operators 
that can be produced, however, wc know surely which operators will be there. 
For simplicity consider unitary minimal models A4{p,p + 1) with the central 
charge c — c and primary operators (prt with conformal weights 

A,, ^ "-'^ i<^<p; i<*<p-i 

When wc have n copies of a CFT wc can label the primary operators with the 
upper indices i like ^^.j; where i = 1,2, n. One can easily show that the OPE 
of two twist fields will not produce one copy of However, (f)l,^(j)l^ which 
is made of two copies of primary operators, with the conformal weight 2 Art, 
always appears after the fusion of two twist operators. In general when we have 
n > 2 the combination of m{l < m < n) copies of the primary operators, with 
the conformal weight mArt, will always appear after the fusion. For n = 2 it 
was shown in |29| that these operators exhaust the fusion structure but for n > 2 
there is no conclusive classification for the fusion structure. In our approximate 
method this will not be a big problem as far as we consider not very big n's 
as we will see soon. The next step is calculating the structure constants. The 
structure constants for small m's were already calculated in and [M]. For 
example for m = 2 one can show 

Clrt = S2(n,4A,,), (28) 



4n2 
i=i 



S2(n,a) = — T2^- (29) 



^This is also related to the modular invariance of the theory. 
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In principle calculating the structure constants for the operators with confor- 
mal weight mA;=rt is related to the calculation of the m-point function of the 
operator ^^t- The calculation is practically doable just for m < 4 and small r 
and < ( to be precise even for to = 4 with r > 2 and t = 1 the calculation is 
very cumbersome). Fortunately, as we will show, most of the terms with m > 2 
give a very small contribution and so one can ignore them in the first approx- 
imation. In addition there are also some cases where the structure constants 
for some fields with to > 2 is zero, this is the consequence of the zero m-point 
function of the corresponding primary operators. For example, since the three 
point function of the spin operators and the energy operators in the Ising model 
is zero, the structure constants for three copies of the spin operator and the 
energy operator are always zero, this is the case also for all odd m's in the Ising 
model. 



1.30 



1.25 



1.20 
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Figure 1: J-2{x) for the Ising model. The blue line is the exact result obtained 
in [17] and the red one comes from the equation (PT|) . 



2.3 Renyi entropy 

In this subsection by using the results of the last two subsections we give an 
approximate formula for the Renyi entropy. Using the formula (|14p with c — nc 
and S ^ ^{1 — ^) one can write the following formula 

nc-l 

^n{x)^ f^^^) " ^Cf(16gr'+..., (30) 

(f3(T))^+--^ V i6g y ^ 

where we take = J-n- Since the above scries is converging very fast, one can 
just take H{c, A;, A^, g) « 1 and also consider just the contribution of the field 
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Figure 2: F'i{x) for the Ising model. The blue line is the exact result obtained 
in [23] and the red one comes from the equation ([51]) . 

with the minimal conformal weight Amm, then 

-.(^) j,;^, (^) * - (5)""^' . ...), (31, 

where we put S2{n) = S2(n, 2A,„i„). Notice that in the above formula is 

two times bigger than the minimum conformal weight in the minimal models 
in (|26l) . To see how good is our approximation we discuss Ising model as the 
simplest unitary minimal model. 

For Ising model and n = 2 the fusion rule is 

T2{zi,zi)T2{z2, Z2)^l + + eh^ ; (32) 

where and cr* are the energy and spin operators of the zth copy respectively. 
In this case Amin = | is the conformal weight of the two copies of the spin 
operator. In Fig. [T]we compare the equation (j3ip with the exact result derived 
in [T7]. The excellent agreement is not a special feature of the critical Ising 
model. Actually we will verify in section |21 where we describe the exact form of 
J^2(a^), that Ea. (|30p for n = 2 is an excellent approximation for any conformal 
model. 

For n = 3 part of the fusion rule is 

73(21,21)73(2:2,22) = 1 + (crV^ + perm) + (e^e^ + perm) + {a^a'^e^ + perm) + (33) 

where ^^perm^' means that we need to consider all the permutations of the three 
indices i = 1,2,3. The dots consider the primary operators that can appear by 
combining the descendants of different copies of the primary operators [31j . 
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Figure 3: The function Fvn{x) for the Ising model as a function of the cross- 
ratio X. The sohd hne is the plot of Eq. while the data are taken form 
Ref. [T7] . The two integers in the legend are respectively the size L of the system 
and the size of the interval i in units of lattice spacing. 



Since the conformal weight of the operator a''(j^e^ is quite bigger than the 
conformal weight of the operator crV-' we do expect that the dominant term 
be still the contribution of &''&■'. In Fig [2] we compare the result coming from 
the equation (|5T|) with the exact result obtained in [21]. Although we ignored 
many contributions, the result is still satisfactory. In principle the result should 
be better for small x's than the large ones, however, for unknown reason our 
formula works better for large x's. It is worth mentioning that if one wants to 
calculate the next terms just considering the terms like a^a^e^ is not enough. 
The reason is that the higher order terms in the expansion of H{c, A;, Ai,q) for 
two copies of a primary field could have equal or bigger contribution. Although 
this is not the case for the Ising model, it simply shows how much it could 
be complicated to go to the higher levels and keep the calculations consistent. 
For future use we also give the first few terms of the result for n > 3. In this 
case the most dominant term after the contribution of the two spin operators is 
the contribution of four spin operators with the conformal weight j. Since the 
power of q in this case in equation (j30]) is one half, one expects an important 
contribution from this term. One can write 
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where S4(n) = Si{n, 2Amin) is calculated in |24| and has the following form 

s^(na) = - V ( sin(7rj42/n) sin(7rj3i/n) V"" ^^^^ 

4 ^ \sin(7rj2i/n)sin(7rj43/n)sin(7rj4i/n)sin(7rj32/n)y 

I<j2<j3<34<"-1 

where jki — jk — ji and ji = 1. Some properties of S2{n^a) and 34(71,0;) are 
discussed in the Appendix. 

2.4 Von Neumann entropy 



Since the equations ([3T|) and (p4| are simple functions of the number of replicas 
n one can simply calculate the von Neumann entanglement entropy by just 
differentiating them with respect to n and evaluating them at n = 1. Then for 
Ising model we have 



4(l)(4g)^+4(l)(4g)^) + ..., ,(36) 

where 33(1) is known from the work in |24) and for generic conformal weight 
has the following form 



(.)(^^j (-6log^^3(.) + ^log(- 



^^['''^™"^^+4r(2A_ + |) ■ 



Calculation of 34(1, -j) is a very difficult task, we find an approximate value for 
this quantity by using a numerical method described in the Appendix which 
yields 5^1(1) ~ 0.126(3). 

We compare Eq. with the numerical data of Ref. [T7] , based on a 
tree tensor network algorithm (TTN) [32] applied to a critical one dimensional 
quantum spin chain in transverse field, corresponding to the CFT minimal model 
with c = i . This algorithm gives the full spectrum of the reduced density matrix 
Pa- From this the moments trp^ and the entanglement entropy 5^1 are easily 
evaluated. The data of [TJ are taken for a periodic system of length L with a 
subsystem A = [vi —ui] U [v2 — U2] composed of two disjoint intervals of identical 
length £ = wi — 1*2 = W2 — 1*2 at distance r. The cross-ratio x is given by 



sin Tr£/L 
sin7r(£ + r)/L 



(38) 



The ratio p\/ pf^ eliminates the non universal constant c„ of Eq. ^ and 
allows to evaluate, up to obvious factors, the universal function !Fn{x) and its 
contribution Fvn = 9„J-"„(a;)|„=i to the von Neumann entropy. In Fig. [3] we 
plot the function J^vn{x) as well as the TTN numerical data, finding a perfect 
agreement. It is interesting to note that if we use the quantity 54(1, j] as a free 
parameter to fit these data we find a value which is consistent, within the errors, 
with the value estimated in the Appendix in a completely different context. 



*We thank the authors of Rcf. 1171 for providing us with their numerical data. 
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3 The special case of two replicas 



The case n = 2 is very interesting and instructive. First, note that Eq. ([50)) 
simplifies dramatically because the structure constants become simply 



l=rt 



4Art 



thus 



:f2{x) 



a;(l — a;) 
16g 



2A, 



(39) 



(40) 



Moreover in this case the four-point function of the twist fields is directly related 
to the formulation of the CFT on a two-sheeted Riemann surface with two cuts. 
This is conformally equivalent to a torus whose modular parameter r is related 
to the position of the branch points by Eq. (|T7l) . As a consequence it is expected 
that the torus partition function Z{t) of the CFT should show up as a factor 

of nziisniisi. 

In this section we would like to reconstruct this relationship within our 
approach by presenting evidence that a factor of Eq. (|40D is just a truncated 
expansion of the partition function Z{t). To this aim, we resort to the following 
useful formulas 



[,(1 _ A = (^.(r)W^4(r))^ 02ir)e,ir)0,ir) = 2,irf . 



(41) 



The Jacobi theta functions 6*2 (r) and 6'3(r) have been already defined in p7|) and 

2 

similarly we have 04(r) — X)nez(~l)"9" i ^i^) ^he Dedekind eta function, 
defined as 



q 



n=l 



These relations, once inserted in Eq. yield 



x{l — x) 
16 



(A,t-c/24)|2 



(42) 



(43) 



where we restored the notation of Eq. and the sum is made over the entries 
of the Kac table of a given minimal model. One recognizes at once that such a 
sum is just the truncation of the first few terms of the partition function 



Z(r) = E |e2'''^('^'-'"'=/24)|2 ^ j^jgj^g^ 



(44) 



thus one is led to conjecture that the exact form of T2{x) is simply 

2c 



x{l — x) 
16 



t]{t) 



3(r) 



Z{r) 



(45) 



Actually this formula coincides for c = 1 with the exact result for compactified 
boson found long time ago P51I551IM] and for c = i with the exact result of the 
critical Ising model obtained in [T7] . 
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4 Conclusion 

In this paper wc found approximate formula (|3ip for the Renyi entropy ( with 
arbitrary n) of two disjoint intervals of generic minimal conformal field theory. 
The idea was based on using the elliptic expansion of the conformal blocks 
appearing in the four point function of the twist fields. Since our formula had a 
simple relation with n, the number of replicas, we were able to do the analytical 
continuation n — >■ 1 easily and found the von Neumann entanglement entropy. 
The von Neumann entropy of two disjoint intervals in the Ising model was 
investigated as a benchmark and we found excellent agreement between our 
formula p6p and the available numerical results. 

Using the perturbative expansion of the four point function of the twist fields 
in the n — 2 case we were also able to find a new way to connect explicitly the 
four point function of the twist fields to the partition function of the conformal 
field theory on the torus. 

Although our perturbative results give rather accurate formulas, even by tak- 
ing just the first term of the expansion, for the Renyi entropy and von Neumann 
entropy, it is still tempting to understand how the results can be improved by 
going to the higher levels. In particular, calculating the von Neumann entropy 
of more complicated minimal models and checking the results with numerical 
calculations can shed some light on the possible generalization of our results. 
Since going to the higher levels of the perturbation is possible just by calculating 
the structure constants appearing in the OPE of two twist fields - which is itself 
related to the k point functions of the primary operators - we are faced with 
the old, unsolved problem of conformal field theory of calculating correlation 
functions for arbitrary number of fields. 
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A The functions S2{n^a) and s^{n^a) 

In this appendix we deal with the functions S2(n, a) and s^(n,a), where a is 
related with the conformal weight by A = a/4. 

The function S2{n^a) has a simple zero at n = 1 while s^in^a) has simple 
zeros at n = 1,2,3. For integer n both functions have moreover a factor n 
produced by the translation invariance on the sum of the indices ji so we may 
assume that both functions have also a zero at n = 0. This is in particular true 
for S2[n, k] for integral fc, where it has been shown in two different contexts and 
two different ways |30j[35] that S2[n,k] is a polynomial of degree 2k. The first 
few polynomials of this set are 




(A.l) 
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(this expression has also been found in 



2 ^ sin''^ V90 ' 9 90 



S2{n, 2)^'-y] ^ „ ( ^ + !L _ 11 1 (A.2) 



, n\X \ n4 An^ 191 \ 

52(71, 3) = - > — 7i — -=71 \ \ , (A. 3) 

2^sin^^ V945 90 45 1890 /' ^ ' 



, 1 / 71^ 471^ 7n4 8 ri^ 2497 \ 

S2(n,4) = -> — s — r=n i \ 1 . (A. 4) 

^ ' ^ 2^sin^^ V945O 2835 675 105 28350/ ^ ' 

It is easy to verify in these cases the general formula of Ref. [53] 



4r(| + a) " 4 \2' 



4(",«)n=l- .^,3 , , =7B(-,1+q), (A.5) 



where B{x,y) is the Euler Beta function. In particular, the contribution of 
the operator energy in critical Ising model is 33(71, 2)„=i = and one could 
add this term to Ea. p6p to further improve the formula for the von Neumann 
entropy. 

Assuming that 52(2, a) and 34(2, a) are analytic functions on the complex 
plane z, they admit the following power expansions 

52(2;, a) = z{z — 1) (co + ciz + C2z'^ + c^z^ + ...), (A. 6) 

S4{z, a) = z{z - l){z - 2)(z - 3) {do + diz + ^2^^ + ^3^^ + ...)• (A. 7) 

If the two series expansions in the last parenthesis are rapidly convergent series 
we have, approximately, 

S2(1,Q!) ~ Co + Ci + C2 + C3, (A. 8) 

and similarly 

4(1, a) ~ 2 (do + di+d2 + da) . (A.9) 

We fitted cq, ci, C2, using the exactly calculable data for 32(11, j) for integer 
n = 2, 3, ... 20 and similarly dQ,di, d2, da; 54(11, i) for integer 11 = 4, 5, ... 20 and 
verified in both cases that the coefficients Ci and di are rapidly decreasing. In 
this way 53(1, j) turns out to be compatible with the exact value. In the same 
way we estimated 

4(1, 1)~ 0.126(3) , (A.IO) 

which is the value we put in Eq. p6p: it turns out to fit nicely the numerical 
data as shown in Fig. [31 Using the fitted parameters Ci and di one could also 
study the analytic continuation of J-"„(j;) to Tz{x), where z is the complex plane 
of replicas [3^ simply by inserting the two functions <\A.6\i and <\A.7\i in Eq. ([3T|) . 
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